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Imaging Atmospheric Cherenkov Telescopes (lACTs) are ground-based instruments devoted to 
the study of very high energy gamma-rays coming from space. The detection technique consists 
of observing images created by the Cherenkov light emitted when gamma rays, or more generally 
cosmic rays, propagate through the atmosphere. While in the case of protons or gamma-rays 
the images present a filled and more or less elongated shape, energetic muons penetrating the 
atmosphere are visualised as characteristic circular rings or arcs. A relatively simple analysis of 
the ring images allows the reconstruction of all the relevant parameters of the detected muons, 
such as the energy, the impact parameter, and the incoming direction, with the final aim to use 
them to calibrate the total optical throughput of the given lACT telescope. We present the results 
of preliminary studies on the use of images created by muons as optical throughput calibrators of 
the single mirror small size telescope prototype SST-IM proposed for the Cherenkov Telescope 
Array. 
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1. Introduction 

The Cherenkov Telescope Array project (CTA) [1] will be the next generation ground based 
observatory in gamma-ray astronomy. It will consist of two arrays of imaging atmospheric Cherenkov 
telescopes (lACTs) installed in the two hemispheres to cover the full sky. In its final configuration, 
the Southern site of CTA will have more than a hundred telescopes of three different sizes arranged 
over an area of the order of few km^. The array has been designed in such a way that it will be able 
to cover the wide energy range between few tens of GeV and at least 300 TeV, with an improved 
sensitivity by about a factor of 10 (at 1 TeV) compared to the existing experiments. The array of 
small-size telescopes (SSTs) will be dedicated to the observation of the most extreme gamma-ray 
sources emitting in energies between a few TeV and 300 TeV. 

The SST-IM project is a Consortium of several Swiss and Polish institutes working on the 
design and construction of one of the prototypes for the SSTs of CTA [2]. The SST-IM pro¬ 
totype is the only SST telescope designed with one mirror. It uses the standard and well-proven 
Davies-Cotton design for the optics and telescope structure and an innovative camera using Silicon- 
Photomultipliers (SiPM), inspired by the experience of the FACT camera [3]. Thanks to their wide 
field-of-view (FoV) of about 9°, the SSTs are particularly well suited for spectral studies of Galac¬ 
tic sources and Galactic Plane Surveys. The prototype of the SST-IM mechanical structure has 
been installed on a test site at the Institute of Nuclear Physics PAS in Krakow in November 2013. 
The design of the camera has been finalized and the commissioning phase has already started [4]. 

The detection technique used by the lACTs consists in observing images created by the Cherenkov 
light emitted by electromagnetic showers induced in the atmosphere by gamma rays, or more gen¬ 
erally cosmic rays. To analyse the observed gamma-ray characteristics, an accurate calibration of 
the camera [6] as well as the evaluation of the instrument optical efficiency need fo be performed. 
Highly energefic muons penefrafing fhe afmosphere are visualised as characferisfic circular rings or 
arcs in fhe camera plane and fhey represenf a powerful and precise mefhod fo calibrafe fhe opfical 
fhroughpuf of lACTs [5]. 

In fhis paper we sfudy fhe possibilify of using fhe analysis of muon ring images as a calibrafor 
for fhe opfical fhroughpuf of fhe SST-IM profofype for fhe Cherenkov Telescope Array. 


2. Muon images in lACTs 

Isolafed muons produce sharply defined ring-like images in fhe focal plane and provide a 
powerful fool for moniforing fhe behaviour of fhe telescope performance characferisfics like fhe 
poinf spread function (PSF) and fhe overall lighf collecfion efficiency. Muons passing fhrough 
fhe cenfre of fhe telescope wifh frajecfories parallel fo ifs opfical axis will produce azimufhally 
symmefric rings in fhe camera. The rings will have radii given by fhe Cherenkov angle of fhe 
muons and fhe fofal number of photons expecfed in fhe ring, mainly relafed fo fhe muon energy, can 
be calculafed from fhe measured value of fhis angle. The correspondence befween fhe Cherenkov 
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angle and the muon energy can be established following the Cherenkov equation: 


cos dc = 


1 


( 2 . 1 ) 


where £"0 ~ 0.105 GeV is the muon rest mass, and n the refractive index of air. The typical ring 
radius (R) is determined by the muon spectrum at the observatory altitude, e.g. 7? ~ 1.2° at 1500 m 
a.s.l. for E^ = 20 GeV. Muons arriving under a certain angle with respect to the telescope’s axis 
will form rings or arcs with centres that are offset from the centre of the camera, while muons 
with non-zero impact parameters' will produce rings or arcs with an azimuthally dependent photon 
density. Since the morphology and location of the muon ring allow the muon’s trajectory and 
energy to be calculated, it is possible to predict with precision the number of Cherenkov photons 
that have been collected by the camera. Since the prediction depends only very slightly on the local 
atmosphere, but goes linearly with mirror reflectivity, transmission of glass coating and the photon 
detection efficiency of the camera pixels, the detector response to local muons can be calibrated in 
absolute terms [7]. 


2.1 Simulations 

The results reported in this contribution have been obtained from simulation studies using 
standard CTA software, namely CORSIKA [8] and the sim_t el array package [9]. The latest 
telescope and camera parameters (e.g. min'or reflectivity, light collectors, photon detection effi¬ 
ciency of the silicon photo-multiplier sensors forming the SST-IM camera, hexagonal geometry of 
the camera, etc.) have been included in the simulation of the instrument. The baseline trigger is a 
digital sum trigger with sampling rate 250 MHz. Muons have been injected at a starting altitude 
of ~ 1 km above the observation level^, randomly distributed in a cone of ~ 4.5° (corresponding 
to half of FoV of the telescope) around the vertical axis to cover the entire FoV. The telescope 
has been simulated pointing always to the zenith, and the impact parameter has been randomly 
distributed from 0 m to ~ 2 m (radius of the mirror). The muon spectrum has been simulated as a 
power law with a spectral index of —2 in an energy range between 6 GeV and 10 TeV An example 
of a full ring image as it would appear on the SST-IM is shown in Figure 1. 

2.2 Analysis of the rings 

The analysis of the muon images requires the calculation of the muon image parameters: the 
ring radius, the ring width (ArcWidth), the opening angle of the ring fraction (ArcPhi) and its total 
light content, measured in photo-electrons (MuonSize). 

The analysis has been carried out using MARS, the standard software for the MAGIC tele¬ 
scopes data analysis [10]. 

The search for muon rings starts after image cleaning. For this study an absolute image clean¬ 
ing has been applied with 5 photo-electrons for core and 2 photo-electrons for boundary pixels [11]. 
After cleaning, every image is fitted by a circle starting at the centre of the Hillas ellipse (Hillas 

*The impact parameter is defined as the distance from the centre of the mirror to the impact point of the muon in 
the plane normal to the optical axis. 

^The reference site for all simulations is Aar, the proposed site for the Southern array in Namibia, located at 
/?tei = 1640 m a.s.l. 
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Figure 1: Example of a muon image recorded by the SST-IM telescope. The azimuthal non-uniformity is 
the result of a non-zero impact parameter, while the offset with respect to the centre indicates an off-axis 
muon direction. 

parameters are calculated in a preceding step). From that starting point, the algorithm calculates 
the distance of the putative image centre to every pixel, and from the distance values the mean 
(weighted by the pixel content) with its deviation (rms) are calculated. The algorithm minimises 
the deviation by changing the coordinates of the assumed centre. Figure 2 shows the resolution of 
the muon ring radius as a function of the energy, defined as {Rrec — Rtrue)/Rtrue- The bias on the 
ring radius reconstruction is less than ~ 1 % for the overall energy range considered and it can be 
easily subtracted. 



Figure 2: Resolution of the muon ring radius as a function of muon energy. The bias on the radius recon¬ 
struction is less than ~ 1%. 

In a second step, the image before cleaning is used to calculate the radial and the azimuthal 
intensity distribution of the image. For the azimuthal intensity distribution, all pixels inside a 
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certain margin around the radius are used. The following parameters, useful for the analysis can be 
defined: 

• ArcWidth is calculated as the sigma value of a Gaussian fit to the signal region in the radial 
intensity distribution; 

• ArcPhi is calculated as the sum of connected bins, which lie above a certain pixel amplitude 
threshold; 

• MuonSize is calculated as the sum of the contents of all pixels along the ring. 

Due to the intrinsic broadening of the ring width, the absolute value of the ArcWidth is not equal 
to the sigma value for the optical PSF, but it is related to that. Therefore, to get the size of the 
PSF it is necessary to compare simulated muons with different PSFs with muons extracted from 
observational data. Figure 3 shows an example of the radial intensity distribution of the image with 
a Gaussian fit to the ring region. 


Width distribution 



Figure 3: Radial intensity distribution with a gaussian fit for the calculation of the ArcWidth. The value of 
the ArcWidth is related to (but not exactly) the optical PSF. 


3. Study of the optical degradation 

The analysis presented in the previous section has been developed to study the capability of 
the SST-IM telescope to trigger on muon events and its ability to calibrate the optical throughput, 
possibly without trigger bias. The effect of the summed pixel amplitude threshold and selection 
cuts have been investigated with the purpose of obtaining a muon efficiency curve independenf 
from fhe opfical degradation of fhe felescope. Simulafions have been carried ouf using differenf 
scales of mirror relleclivily degradation. In order fo selecf good muon evenfs^ fhe following qualify 
cufs have been applied based on fhe previously described analysis: 

^Cuts are based on muon simulations only. No study about the cosmic ray rejection has been performed at this 
stage. 
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• the ring fitting the muon image is fully contained in the camera; 

• the reconstructed radius is between 0.5° and 1.5°; 

• the projected ring width along the ring radius has a gaussian shape (reduced chi square of the 
ArcWidth fit < 2); 

Figure 4 (left) shows the muon efficiency (calculated as the ratio between the observed and simu¬ 
lated muons) curves obtained with a trigger threshold of 100 ADC counts as a function of the muon 
energy above 10 GeV. The effect of the muon energy on the muon efficiency is rather negligible 
for all the considered optical efficiencies. The right side of Figure 4 shows the muon efficiency 
for energies above 10 GeV as a function of the optical efficiency. The muon efficiency has a con¬ 
stant value of around 60% in the region of interest (between 60% to 100% of optical efficiency). 
The trigger rate seems to be rather insensitive to the optical efficiency of the system when a lower 
trigger threshold of 100 ADC counts is adopted^. 
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Figure 4: Left plot: Muon efficiency for SST-IM obtained with a trigger threshold of 100 ADC counts 
as a function of energy for different optical efficiencies (100% in red, 90% in green, 80% in blue, 70% in 
yellow, 60% in magenta, and 50% in cyan). Right plot: Muon efficiency for SST-IM obtained with a trigger 
threshold of 100 ADC counts and as a function of the optical efficiency for muons with an energy above 10 
GeV. 


4. Conclusions 

Muon rings have been used as a method to calibrate the total optical throughput of practically 
all previous Imaging Atmospheric Cherenkov Telescopes [12, 13, 14]. In this contribution we 
have presented simulations and analysis of the ring images carried out in the framework of the 
SST-IM project. This work is done in close collaboration with the other SST teams of the CTA 
Consortium [15, 16]. Preliminary results show that muon images are adequately selected despite 
a degradation of the optical efficiency of the instrument. Building on this encouraging result the 
SST-IM collaboration continues to work towards further improving the trigger efficiency for muon 
rings. In this direction, additional dedicated muon triggers are currently under study. 

^the general Prod2 simulation for CTA has been done using a trigger threshold of 127 ADC counts. 
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